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SUMMARY 


INTRODUCTION 

The  development  and  performance  evaluation  of  ASW  sonobuoys  requires  analysis  of  the 
sonobuoy  interaction  with  the  ocean  environment.  Ocean  currents  and  waves  affect  the  sensor 
depth,  orientation,  self-noise,  and  survivability  and  must  be  taken  into  account  in  the  earliest  phases 
of  sonobuoy  design  and  development  An  accurate  computer  simulation  permits  early  analysis  of 
proposed  systems,  provides  a  basis  for  evaluating  changes  in  the  system  as  the  design  evolves,  and 
establishes  a  tool  for  analyzing  performance  of  component  and  prototype  hardware  during  test  and 
evaluation. 

Perhaps  the  most  useful  computer  program  developed  for  analysis  of  sonobuoy  systems  is 
entitled  CABUOY.  This  program  has  the  capability  of  modelling  a  frae-floating,  moored,  or  towed 
system  in  a  two  dimensional  ocean  current  with  wave  excitations  from  sinusoids  to  random  sea 
states.  The  complexity  of  a  computer  program  like  CABUOY  requires  an  extensive  validation  to  be 
undertaken.  It  is  only  by  performing  experiments  and  comparing  the  computer  results  with  an 
empirical  “ground  truth"  that  confidence  in  the  accuracy  of  the  computer  program  can  be  ob¬ 
tained.  To  this  end,  a  set  of  validation  experiments  were  conducted  to  verify  computer  results. 

This  report  presents  the  results  of  the  first  experimental  validation  of  the  CABUOY  computer 
program  in  the  simplest  case  considered,  i.e.,  the  vertical  motion  of  a  body  suspended  in  still  water 
at  the  lower  end  of  an  elastic  cable  when  the  upper  end  of  the  cable  (the  surface  float  end)  was  sub¬ 
jected  to  vertical  harmonic  motion. 

BACKGROUND 

The  dynamic  representation  in  a  computer  program  for  an  ocean  deployed  system  consisting 
of  a  surface  buoy,  connecting  cables,  and  one  or  more  submerged  body  has  been  considerable 
interest  to  the  oceanographic,  sonobuoy,  and  sonar  communities.  A  time  domain  computer  pro¬ 
gram  entitled  CABUOY  was  developed  for  the  Naval  Air  Development  Center  (NAVAIRDEVCEN) 
by  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC)  to  model  this 
ocean  deployed  buoy  system  1/2,3. 

Because  of  CABUOY's  utility,  immediate  use  was  made  of  the  program  to  model  sonobuoy 
systems.  As  a  result,  several  modifications  were  implemented  by  the  NAVAIRDEVCEN  and  the 
DTNSRDC  to  expand  or  improve  the  program4,5,6,7. 

RESULTS 

The  validation  of  the  CABUOY  program  for  the  vertical  motion  of  a  spring-body  suspension 
system  in  water  was  successful  but  only  after  some  revision  was  made  to  the  program  to  correct  an 
iteration  calculation  of  virtual  mass. 

The  internal  calculations  of  the  (revised)  CABUOY  program  for  virtual  mass  and  drag  of  disks 
were  required  for  valid  results  of  systems  with  disks.  A  single  disk  and  two  independent  disks  in  the 
same  suspension  were  successfully  modelled  by  CABUOY  using  the  internal  routine  for  virtual  mass. 

The  CABUOY  modelling  of  a  system  with  a  low-aspect  ratio  cylinder  (and  no  disk)  was  only 
approximate  because  of  the  approximation  of  virtual  mass  and  drag  for  this  body. 
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CONCLUSIONS 

The  accuracy  of  the  CABUOY  program  results  depends  upon  correct  user  Input  to  characterize 
the  spring,  especially  a  non-linear  spring  normally  used  in  a  sonobuoy  system,  where  attention  must 
be  paid  to  the  dynamic  as  well  as  static  spring  properties. 

Further  data  on  virtual  mass  of  oscillating  bodies  is  required  to  permit  the  user  a  wider  scope 
of  CABUOY  validity. 

RECOMMENDATIONS 

The  body  of  data  on  the  virtual  mass  and  drag  of  oscillating  bodies  of  various  geometries  needs 
to  be  expand^  for  across-the-board  utilization  of  the  CABUOY  program. 

A  modification  to  the  CABUOY  program  to  permit  a  more  direct  input  of  spring  character¬ 
istics  should  be  made.  This  modification  might  employ  a  separate  static  and  dynamic  spring 
constant 

A  revised  program  listing  and  additionai  user's  infomnation  should  be  provided  to  all  those 
users  who  have  the  CABUOY  program. 

Other  aspects  of  CABUOY  should  be  validated. 
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DISCUSSION 


CABUOY  SYSTEM  MODEL 

CABUOY  solves  for  the  dynamic  motions  of  drifting  (free  floating),  moored,  or  towed  systems 
consisting  of  a  surface  float  and  any  combination  of  suspended  bodies  and  interconnecting  cables  as 
shown  in  figure  1.  The  surface  float  can  be  input  as  a  spheroid,  ranging  from  a  thin  disk  (extreme 
oblate  spheroid)  to  a  spar  (extreme  prolate  spheroid),  which  interacts  with  the  surface  waves,  or  the 
float  can  be  directly  coupled  to  or  completely  decoupled  from  the  surface  for  prescribed  motions 
under  experimental  conditions.  The  cable  is  divided  into  a  number  of  rigid  (straight  line)  extensible 
segments.  A  body  represented  by  a  point  mass  of  known  drag  area  and  virtual  mass  may  be  in¬ 
cluded  between  each  cable  segment,  if  desired. 

The  ocean  environment  includes  a  user-defined  sea  surface  which  decays  exponentially  with 
depth.  This  is  generated  by  the  summation  of  up  to  nineteen  sine  waves  with  various  amplitudes, 
frequencies,  and  phases.  The  program  will  also  internally  gen;<rate  sea  states^.  A  two  dimensional 
current  profile  of  coplanar  horizontal  current  velocity  with  water  depth  is  input,  as  well.  The 
steady-state  calculations  for  the  catenary  in  the  current  profile  use  the  finite  element  analysis 
techniques  and  interation  routines  used  in  the  steady-state  free  floating  and  moored  system  com¬ 
puter  programs  FF2E  and  MR3E8.  9. 

The  cable  stress-strain  relationship  is  defined  by: 

'*’  =  "I'ref*‘Cl  €^2  +  C|(^je  eq(1) 


where 

T  =  cable  tension 

Tref  =  constant  "reference"  tension 

Cl,  C2  =  constants 

C|NT  ~  internal  damping  coefficient 

€  =  Strain,  i.e.,  change  in  length/original  length 

i  =  rate  of  change  of  strain,  de/dt 

For  most  applications  Cinj  is  small  and  the  last  term  in  equation  (1)  can  be  neglected.  If  the  cable 
is  a  linear  spring,  C2  =  1  and  Ci  >  AE,  where  A  is  the  cable  cross  sectional  area  and  E  is  the  elastic 
modulus.  A  linear  approximation  of  the  cable  is  generally  made  around  the  operating  tension,  as 
illustrated  in  figure  2.  Tref  can  be  found  from  the  T-axis  intercept  of  the  linearized  curve  and 
Cl  *  dT/de,  or  the  slope  of  the  linear  portion  of  the  curve.  The  computer  program  uses  this  equa¬ 
tion  to  determine  the  static  elongation  of  the  cable  to  relate  the  dynamic  change  of  tension  to  the 
amplitude  of  the  motion  of  the  body  (or  bodies)  attached  to  the  C8ble(s). 


k 


,V> 


3 


NAOC-83112-30 


The  hydrodynamic  effects  on  a  body  being  moved  through  a  fluid  are  two-fold.  A  drag  force 
which  resists  the  motion  is  velocity  dependent  as  described  by  the  equation: 

Fd=7CdSV2  eq(2) 

where 

Fq  *  drag  force 
p  *  fluid  density 
Cq  s  drag  coefficient 
S  »  cross-sectional  area  of  body 

V  «  velocity  of  the  body  relative  to  the  fluid 

If  the  body  is  moving  at  constant  velocity,  the  drag  force  (assuming  no  lift  force  applies)  is  the  only 
hydrodynamic  effect  of  concern.  However,  if  the  body  is  accelerated  through  the  fluid,  a  second 
effect  is  of  importance,  i.e.,  the  hydrodynamic  mass.  An  accelerating  body  has  an  effective  mass, 
called  the  virtual  mass,  which  is  defined  as 

mv=mp  =  mh  eq  (3) 

where 

my  =  virtual  mass 

nnp  =  physical  mass  =  weight  (in  air)/gravitational  acceleration 
mh  =  hydrodynamic  (or  added)  mass. 

The  hydrodynamic  mass  is  more  elusive  in  character  than  drag  because  it  is  more  difficult  to  mea¬ 
sure,  but  it  becomes  a  significant  factor  in  determining  the  motion  of  a  body  oscillating  in  a  fluid. 
Included  as  an  option  in  the  CABUOY  program  is  an  internal  calculation  of  the  hydrodynamic  mass 
for  disks.  Based  on  a  previous  empirical  investigation^^,  this  calculation  was  an  iterative  operation, 
determining  the  hydrodynamic  mass  and  drag  of  a  disk  in  harmonic  acceleration  from  the  peak 
values  of  velocity  and  acceleration  according  to  the  following  equations: 

mh  ■  1.2pd3  (v2/ad)^/2  eq  (4) 

Co  -  2.2  (v2/ad)-V2  eq  (5) 

where 

mn  ■  hydrodynamic  mass 

p  >  fluid  density 

d  >  disk  diameter 

V  ■  peak  harmonic  velocity 
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a  -  peak  harmonic  acceleration 
Cq  *  drag  coefficient 
This  is  shown  in  figures  3  and  4. 

The  CABUOY  program  prints  out  the  horizontal  and  vertical  positions  of  the  surface  float  and 
each  cable  node  or  body  referenced  to  the  surface  float,  their  vertical  and  horizontal  velocity  and 
acceleration,  their  tilt  angle  from  vertical,  and  the  tension  in  the  cable  at  the  body  for  each  point  in 
time,  spaced  at  user-selected  time  intervals. 

VERTICAL  MOTION  VALIDATION  MODEL 

The  primary  experiment  to  evaluate  the  CABUOY  vertical  motion  dynamics  employed  a  body 
suspended  from  an  elastic  cable.  The  system  was  »jbmerged  in  still  fresh  water  in  the  Naval  Su'  e 
Weapons  Center  (NAVSURFWPNCEN)  100-foot  deep  Undersea  Weapons  Tank,  as  shown  in 
figure  5.  A  rotating  arm  drive  was  used  to  produce  a  harmonic  motion  forcing  function  which 
translated  to  vertical  motion  by  means  of  a  pulley.  Input  amplitude  and  frequency  could  be  va  . 
at  the  motor.  The  motion  of  the  body  was  view^  through  portholes  against  a  background  grir  I 
recorded  by  an  observer  as  well  as  photographed  by  a  motion  picture  camera.  In  addition,  a  Ic 
frequency  response  pressure  transducer  system  calibrated  to  give  vertical  amplitude  of  the  bod  ■' 
recorded. 

The  bodies  used  were  of  three  types,  as  shown  in  figure  6.  A  3-inch  high,  4-inch  diameter 
aluminum  cylinder  weighing  4.24  lbs  in  air  and  2.92  lbs  in  water  was  the  basic  body.  To  this,  was, 
at  times,  att^ed  a  12-inch  diameter  disk  of  1/4'inch  thick  plexiglas  (0.96  lbs  in  air,  0.16  lbs  in 
water),  a  12'inch  diameter  plastic  flexible  disk  (of  the  type  used  in  sonobuoy  dampers,  e.g.,  a  spring 
steel  hoop  stretching  taut  a  few  mil  thick  plastic  sheeting),  and  a  pair  of  12-inch  plexiglas  disks 
separated  by  5-feet  of  inextensible  line. 

The  elastic  cable,  or  spring  of  this  spring-mass  system,  was  a  30-foot  long  piece  of  multistrand 
bungee  of  the  type  us^  in  several  sonobuoy  systems.  Typically,  the  bungee  is  a  non-linear  springl  1 
and  measurements  of  the  bungee  were  made  to  characterize  it  for  input  to  the  computer  program. 
The  details  of  the  bungee  characterization  will  be  discussed  later.  Initially,  a  simple  stretching  of 
the  bungee  was  used  to  generate  a  tension-elongation  plot  as  shown  in  figure  7. 

SINGLE  DISK  SYSTEM 

The  single  disk  system  (figure  6b)  was  selected  as  the  first  subject  of  comparison  between  mea¬ 
sured  behavior  and  computer  predicted  performance  because  the  internal  calculation  of  virtual  mass 
could  be  utilized.  The  experimental  data  are  shown  in  table  I.  Runs  6-12  used  a  rigid  plexiglas  disk 
and  Runs  13-19  used  a  flexible  plastic  disk.  No  significant  difference  between  the  two  disks  were 
noted.  Initial  CABUOY  outputs  bore  little  resemblance  to  the  experimental  results.  In  an  effort  to 
determine  why  the  computer  results  differed  from  empirical  data,  the  program  was  modified  to 
print  out  the  body  virtual  mass  and  drag  coefficient  at  each  time  interval  in  the  output.  The  values 
that  the  progrwn  was  calculating  were  not  changing;  they  should  have  with  the  iterative  process  the 
CABUOY  program  was  supposed  to  perform.  The  CABUOY  generated  values  corresponded  to 
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those  expected  at  the  limits  placed  upon  the  hydromechanical  mass  and  drag  coefficients  of  equa¬ 
tions  (4)  and  (5),  i.e., 

0.08  <v2/ad<  3.0.  eq  (6) 

The  parameters  V  and  a  were  apparently  attaining  values  which  would  cause  V^/ad  to  achieve  one 
of  the  limits. 

The  lower  limit  on  v2/ad  is  based  on  small  amplitude  vibration.  In  harmonic  oscillation  the 
peak  velocity  is  the  peak  amplitude  of  motion,  y,  multiplied  by  the  angular  frequency,  co: 


< 

It 

< 

E 

eq  (7) 

and  similarly 

a  *  y6j2 

eq  (8) 

therefore. 

v2/ad  =  y2Gj2/y^2c|  =  y/<j 

eq  (9) 

The  acceleration  modulus  V^/ad  is  a  measure  of  the  amplitude  of  harmonic  motion  in  disk  dia¬ 
meters.  At  small  amplitudes,  the  disk  can  be  considered  a  vibrating  circular  piston  in  an  infinite 
battle12  which  for  small  amplitudes  yields  a  hydrodynamic  mass  approaching 

f'^h-jpcl^"  eq  (10) 

This  is  the  same  value  as  the  theoretical  hydrodynamic  mass  from  the  kinetic  energy  imparted  bv  a 
moving  disk  to  an  ideal,  incompressible  fluid  of  infinite  extent  as  calculated  by  Milne-Thompsonl3. 

As  the  amplitude  of  motion  increases,  a  vortex  ring  which  forms  on  the  disk  influences  the 
hydrodynamic  mass.  At  large  amplitudes,  the  vortex  sheds  and  the  disk  can  be  caused  to  skew  to 
the  side,  "dumping"  entrained  water  and  confusing  the  dynamics  to  the  point  where  repeatability 
is  lost  This  is  the  reason  for  the  upper  limit  of  v2/ad.  Even  with  these  bounds,  a  12-inch  dia¬ 
meter  disk  should  have  its  hydrodynamic  mass  defined  by  CABUOY  for  a  peak-to-peak  motion 
amplitude  between  1.92  inches  and  72.0  inches. 

Removing  the  limits  in  the  program  did  not  resolve  the  problem.  Further  study  of  the 
CABUOY  program  uncovered  a  misinterpretation  of  V  and  a  in  the  virtual  mass  calculation.  In¬ 
stantaneous  valves  were  being  used  instead  of  peak  values.  CABUOY  was  modified  to  retain  maxi¬ 
mum  valves  for  V  and  a  over  a  time  interval  and  use  these  values  ir  recomputing  the  motion.  This 
would  be  repeated  during  the  next  computation  interval,  and  so  on.  This  iterative  technique  finally 
resulted  in  a  steady-state  amplitude  with  consistent  virtual  mass  calculations. 

SPRING  CHARACTERISTICS 

The  compliant  cable  used  in  most  sonobuoy  suspensions  to  isolate  the  sensor  from  the  ocean 
wave  effects  on  the  surface  float  is  a  synthetic  or  natural  rubber  elastomer  (or  bungee  cord).  The 
bungee  cord  used  in  the  validation  testing  was  multistranded  (9-strand)  rubber  with  interwoven 
thread  to  provide  a  "haired-fairing"  to  reduce  strumming  in  relative  flow. 
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According  to  Hookes'  Law, 


T=  ksA6=  ks  (C-fio) 


eq{11) 


where  T  is  the  force  (tension)  applied,  S  is  the  stretched  spring  length  under  tension  T,  £o  is  the  un¬ 
stretched  spring  length,  and  ks  is  the  (static)  spring  constant.  The  spring  constant  can  t)e  found  by 
tensioning  the  spring  by  increasing  amounts  and  measuring  the  resultant  stretch.  By  setting 
elongation 


eq(12) 


T  -  ( ks£Q)  € 


eq  (13) 


which  is  the  form  of  equation  (1)  where  ks£o  =  Cf,  Tref C|NT  =  0,  and  C2  =  1.  Because  a  bungee 
cord  is  not  a  linear  spring  but  exhibits  behavior  like  that  shown  in  figure  2,  a  portion  can  be  linear¬ 
ized  around  the  operating  point  by  using  equation  (1).  This  is  shown  graphically  in  figure  7  for 
the  multistrand  bungee  used. 


Another  way  of  determining  a  spring  constant  is  through  the  solution  to  the  equation  of 
motion  of  an  undamped  spring-mass  system: 


my  ky  *  0 

The  solution  to  this  equation  is 


eq(14) 


y  =  yoSincAjt  eq  (15) 

where  yo  is  the  peak  amplitude  of  motion,  t  is  time,  and  u  is  the  angular  frequency.  The  natural 
frequency  is 


a»n=NrR57m  eq  (16) 

where  kd  is  the  (dynamic)  spring  constant  and  m  is  the  mass  of  body  on  the  end  of  the  spring.  In 
an  ideal  spring, 

ks=kd;  eq  (17) 

however,  this  is  not  true  for  the  bungee  cord.  When  masses  of  different  weight  were  hung  upon  a 
sample  of  the  bungee  and  the  natural  frequency  measured,  kd  was  found  for  the  range  of  loads  on 
the  cable.  This  is  shown  in  figure  8  in  comparison  with  ks  measured  from  the  static  stretch  of  the 
bungee. 

The  linear  static  equation  (as  shown  in  figure  7)  which  satisfies  the  form  of  the  CABUOY  pro¬ 
gram  (eq  (1))  in  the  operating  range  of  loads  (static  load  2.9  to  3.1  lbs)  is 

T -1.08+ 1.05  e  eq(18) 
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In  obtaining  this  equation,  a  graphical  method  was  used  as  illustrated  in  figure  2.  A  similar  equation 
for  the  dynamic  properties  of  the  bungee  can  be  devised.  The  derivation  of  that  equation  is  ex¬ 
plained  as  follows: 

The  static  spring  constant  kg  is  defined  as 


where  Ts  is  the  static  tension,  £o  the  unstretched  length  of  the  spring,  and  £$  is  the  stretched 
length  under  static  tension,  Ts.  This  can  be  written  as 


®oks  “  Ts/cj 

eq(20) 

where 

£s  —  fio 
£o 

eq(21) 

The  non-linear  bungee  has  a  linear  range  described  by  the  equation 

Ts  =  Tref  ^s- 

eq(22) 

Combining  this  with  equation  (20), 

£oks  «s®  "Cref  Cl  «$ 

eq(23) 

Cl  =  £oks  -  Tref/e, 

(eq  (24) 

The  dynamic  spring  constant  can  be  defined^^  as 


A8  £  -  fi. 


eq  (25) 


where  T  is  the  tension  and  £  is  the  length  of  the  spring  displaced  from  the  statis  position.  The 
tension  can  be  expressed  as 


T  =  Ts+kd(£-£s). 


eq(26) 


Since 


£-£o-{fi-£,)  +  (£s-£o).  eq(27) 

equation  (26)  can  be  rewritten  as 

T-T$  +  kd(£-£o)-kd(*s-8o)-  (28) 

Using  equation  (22)  to  substitute  for  Ts, 

T  -  Tref  +  Ci  e,  +  kd  (£  -  £o)  -  kd  («s  -  «o>-  <29) 
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Substituting  equation  (24), 

T  *  T|igf  +  Sgks  “  "^ref  ~ 

or 

T  —  Eq^s  +  kd  ~  ®o^  ~  ^  ~  ^o^- 

Using  equation  (21)  and  rearranging, 

T  =  ks  (Bs  -  fio)  -  kd  {«s  -«o)  +  kd  («  -  «o) 

T  -  (ks  -  kd)  (fis  -  «o>  +  kd  («  -  80)- 
Multiplying  by  Bq^^o* 

fi  —  B 

T  *  (kj  —  kd)  (B5  —  Bq)  +  Bq  kd  «  ^ 

*0 

This  is  of  the  form  of 

T  =  Tref  +  C)  € 

where 


Cl  *  Bq  kd 


and 

Tfgf  •  (kj  —  kd)  (Bj  —  Bq) 
which  by  using  equation  (19)  reduces  to 


eq(30) 

eq  (31) 

eq  (32) 
eq(33) 

eq(34) 

eq  (35) 

eq(36) 

eq(37) 

eq(38) 


Tref“<ks-  kd) 


Jj. 

k| 


eq(39) 


The  dynamic  equation  to  be  used  in  CABUOY  is  then  equation  (35)  with  the  constants  and 
parameters  defined  in  equations  (36),  (37),  and  (39).  For  the  multistrand  bungee  where  Bq  «  30  ft, 
at  Ts  «  2.9  lbs,  ks  ^  0.0567  and  kd  ae  0.0425  and  equation  (35)  becomes 


^  T -0.75+  1.27  e  eq  (40) 

A  comparison  between  the  dynamic  and  static  expressions  describing  the  cable  is  shown  in  figure  9. 
Both  lines  go  through  the  static  operating  point,  but  have  different  slopes.  Both  equation  ( 18)  and 
equation  (40)  were  used  in  CABUOY  to  describe  the  bungee  and  the  results  compared  with  ob¬ 
served  data,  as  shown  in  table  II  and  figure  10.  The  dynamic  spring  equation  produced  slightly 
better  results  than  the  static  equation.  The  CABUOY  results  using  the  dynamic  spring  character¬ 
ization  were  consistent  with  consistently  lower  amplitudes  than  the  observed  values.  A  maximum 
error  of  20%  was  observed  but  most  calculated  data  were  within  12%  of  the  observed  results  except 
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at  small  amplitudes  where  the  largest  experimental  percentage  errors  occurred.  In  absolute  terms, 
the  error  amounted  to  approximately  4  inches  in  34  inches,  2.5  inches  in  13  inches,  and  0.5  inches 
in  1.5  inches,  while  the  visual  observations  had  probable  errors  of  ±1.5  inches. 

AMPLITUDE  VARIATION 

A  phenomenon  observed  in  the  past  was  noted  in  the  vertical  motion  experiments.  When  com¬ 
paring  the  magnification  factor  (i.e.,  output  amplitude/input  amplitude)  of  the  same  system  with 
two  different  input  amplitudes,  the  magnification  factor  of  the  smaller  input  was  larger.  The  disk 
system  with  a  3-foot  input  amplitude  (table  Hi)  was  compared  with  the  6-foot  input  data  (table  I) 
with  similar  results.  Table  IV  shows  the  CABUOY  calculated  amplitudes  for  the  3-foot  input,  and 
figure  1 1  shows  the  comparison  of  observed  and  calculated  data  for  the  3-foot  and  6-foot  input 
cases.  The  magnification  factor  for  the  CABUOY  calculations  follows  the  form  of  the  observed 
data  closely. 

Figure  12  compares  the  CABUOY  calculated  data  for  both  amplitude  inputs  to  the  disk  system 
to  observed  data.  If  the  maximum  error  in  observation  is  ±1.5  inches,  the  calculated  data  falls  with¬ 
in  the  envelope  of  expected  values  only  if  an  offset  in  the  calculation  is  taken  into  account  The 
calculated  values  are  apparently  10%  lower  dian  the  experimental  data.  This  is  probably  a  function 
of  the  coefficients  used  in  the  virtual  mass  and  drag  equations  for  the  disk.  An  area  for  further 
investigation  might  be  in  varying  these  coefficients  for  improved  agreement  with  the  observed  data. 

TWO  DISK  SYSTEM 

The  two-disk  system  (figure  6c)  was  oscillated  with  a  3-foot  input  amplitude  vertical  motion 
on  the  30-foot  bungee.  The  disks  were  separated  by  5  feet  of  non-compliant  line.  Observed  and 
calculated  data  are  shown  in  table  V.  Each  disk  was  considered  as  a  separate  body  in  the  CABUOY 
simulation  with  independently  calculated  virtual  masses  (i.e.,  assuming  no  interference  effects),  as 
predicted  by  reference  10  for  this  separation. 

The  results  show  excellent  agreement  witii  observations,  better  than  the  single  disk  system.  A 
comparison  between  the  single  and  dual  disk  systems  with  the  same  input  amplitude  is  shown  in 
figure  13. 

CYLINDER  SYSTEM 

The  3-inch  high,  4-inch  diameter  cylinder  in  figure  6a  was  oscillated  on  the  end  of  the  30-foot 
bungee  with  both  3-foot  and  6-foot  peak-to-peak  input  amplitudes.  The  CABUOY  program  can 
internally  calculate  the  virtual  mass  and  drag  associated  with  an  oscillating  disk  in  the  form  of 
equations  (4)  and  (5).  Figures  14  and  15  show  the  drag  and  virtual  mass  for  oscillating  cylinders 
over  a  limited  range  of  length-to-diameter  ratios. '*0  Since  the  cylinder  used  in  the  experiments 
had  an  S/d  ratio  of  0.75,  the  disk-type  drag  (see  figure  14)  could  be  used,  but  a  modification  to  the 
virtual  mass  equation  u^  in  the  CABUOY  program  was  required,  since  the  internal  virtual  mass 
calculations  were  originally  only  defined  for  disks. 

The  hydrodynamic  mass  for  the  cylinder  in  figure  15  is  in  the  form 

mh-C|p.Jd28  ap(41) 
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Rather  than  to  introduce  a  new  set  of  equations  into  CABUO  Y,  the  form  of  equation  (4)  was  used; 


mh  =  C|  p^d2  S  =  C*  p  d3  (v2/ad)^/2 

eq(42) 

or 

C|  =  f  fe)c|  (v2/ad)-1/2 

eq(43) 

For  E/d  *  0.75, 

C*  =  0.589  C|  (v2/ad)-^/2 

eq  (44) 

Using  figure  15,  values  for  C|  and  V^/ad  were  found  to  yield  an  approximate  solution  C*{  =  0.55, 
for  £/d  -  0.75.  This  was  used  in  place  of  the  1.2  value  in  equation  (4)  and  the  output  data  com¬ 
pared  with  observed  data  in  figure  16  and  table  VI.  The  agreement  between  the  calculated  and  ob¬ 
served  data  is  not  as  good  as  that  for  the  disk  system  but  is  certainly  a  useful  approximation. 

It  was  beyond  the  scope  of  the  present  study  to  extend  the  virtual  mass  data  for  oscillating 
cylinders,  but  it  became  obvious  that  the  validity  of  the  computer  solutions  are  dependent  on  this 
data.  A  preliminary  graph  of  the  hydrodynamic  mass  coefficient  Cf  in  the  equation 

mhC*  pd3(v2/ad)1/2  eq  (45) 

for  low  S/d  cylinders  is  shown  in  figure  17,  based  on  the  curves  in  figure  15.  An  approximation  for 
drag  coefficient  for  0.8  <  S/d  <  2  may  be  used  as  follows: 

Co  -  0.83  +  0.25  lv2/ad)-1/2  eq  (46) 

whereas  for  S/d  <  0.75,  the  disk  formulation  can  be  used  for  drag. 
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TABLE  I  -  EXPERIMENTAL  DATA  FOR  12-INCH  DIAMETER  DISK  ON  30-FEET  LENGTH 
OF  BUNGEE.  (RUNS  6-12  FOR  PLEXIGLAS  (RIGID)  DISK,  RUNS  13-19  FOR  PLASTIC 
(FLEXIBLE)  DISK)  WITH  6-FEET  INPUT  AMPLITUDE. 


RUN 

PERIOD 

(sec) 

FREQUENCY 

(Hz) 

INPUT 

AMPLITUDE 

PEAK-TO-PEAK 

(ft) 

OUTPUT 

AMPLITUDE 

PEAK-TO-PEAK 

(ft) 

MAGNIFICATION 

FACTOR 

yout/yin 

6 

3.6 

0.278 

6.00 

0.125 

0.021 

7 

5.0 

0.200 

6.00 

0.250 

0.042 

8 

8.0 

0.125 

6.00 

0.474 

0.079 

9 

12.3 

0.081 

6.00 

0.639 

0.107 

10 

14.9 

0.067 

6.00 

1.062 

0.177 

11 

19.6 

0.051 

6.00 

1.437 

0.240 

12 

30.3 

0.033 

6.00 

2.811 

0.469 

13 

3.6 

0.278 

6.00 

0.114 

0.019 

14 

5.3 

0.188 

6.00 

0.182 

0.030 

15 

8.0 

0.125 

6.00 

0.400 

0.067 

16 

12.0 

0.083 

6.00 

0.692 

0.115 

17 

15.9 

0.063 

6.00 

1.120 

0.187 

18 

22.7 

0.044 

6.00 

1.832 

0.305 

19 

30.3 

0.033 

6.00 

2.744 

0.457 
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TABLE  II  -  COMPARISON  OF  CABUOY  CALCULATED  OUTPUT  USING  STATIC  AND 
DYNAMIC  SPRING  EQUATIONS  WITH  OBSERVED  AMPLITUDE  MEASUREMENTS  FOR 
THE  DISK  SYSTEM  (6-FEET  INPUT  AMPLITUDE). 

STATIC  SPRING  DYNAMIC  SPRING 

EQUATION  EQUATION 

OBSERVED  CALCULATED  CALCULATED 


RUN 

FREQUENCY 

(Hz) 

OUTPUT 

(ft) 

OUTPUT 

(ft) 

OUTPUT 

(ft) 

6 

0.278 

0.125 

0.09 

0.08 

7 

0.200 

0.250 

0.16 

0.16 

8 

0.125 

0.474 

0.34 

0.34 

9 

0.081 

0.639 

0.63 

0.63 

10 

0.067 

1.062 

0.83 

0.85 

11 

0.051 

1.437 

1.24 

1.27 

12 

0.033 

2.811 

2.39 

2.49 

13 

0.278 

0.114 

0.10 

0.10 

14 

0.188 

0.182 

0.19 

0.20 

15 

0.125 

0.400 

0.34 

0.34 

16 

0.083 

0.692 

0.61 

0.63 

17 

0.063 

1.120 

0.92 

0.93 

18 

0.044 

1.832 

1.56 

1.60 

19 

0.033 

2.744 

2.41 

2.51 

TABLE  III  -  EXPERIMENTAL  DATA  FOR  12-INCH  DIAMETER  DISK  ON  30-FEET 
BUNGEE  WITH  3-FEET  INPUT  AMPLITUDE. 


RUN 

PERIOD 

(sec) 

FREQUENCY 

(Hz) 

INPUT 

AMPLITUDE 

P-P 

(ft) 

INPUT 

AMPLITUDE 

P-P 

(ft) 

MAGNIFICATION 

FACTOR 

yout/yin 

27 

3.5 

0.286 

3.00 

0.075 

0.025 

28 

6.0 

0.167 

3.00 

0.140 

0.047 

29 

10.0 

0.100 

3.00 

0.321 

0.107 

30 

15.0 

0.067 

3.00 

0.576 

0.192 

31 

20.0 

0.050 

3.00 

0.885 

0.295 

32 

30.0 

0.033 

3.00 

1.749 

0.583 
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TABLE  IV  -  COMPARISON  OF  CABUOY  CALCULATED  OUTPUT  (USING  DYNAMIC 
SPRING  EOUATION)  WITH  OBSERVED  MEASUREMENTS  FOR  DISK  SYSTEM  (3-FEET 

INPUT  AMPLITUDE). 


RUN 

FREQUENCY 

(Hz) 

OBSERVED 

OUTPUT 

P-P 

(ft) 

CALCULATED 

OUTPUT 

P-P 

(ft) 

27 

0.286 

0.075 

0.050 

28 

0.167 

0.140 

0.140 

29 

0.100 

0.321 

0.300 

30 

0.067 

0.0576 

0.520 

31 

0.050 

0.885 

0.820 

32 

0.033 

1.749 

1.640 

TABLE  V  -  COMPAIRSON  OF  CABUOY  CALCULATED  OUTPUT  WITH  OBSERVED 
MEASUREMENTS  FOR  A  TWO-DISK  SYSTEM. 


RUN 

PERIOD 

(sec) 

FREQUENCY 

(Hz) 

INPUT 

AMPLITUDE 

P-P 

(ft) 

OBSERVED 

OUTPUT 

AMPLITUDE 

P-P 

(ft) 

CALCULATED 

OUTPUT 

AMPLITUDE 

P-P 

(ft) 

33 

4.0 

0.250 

3.00 

0.063 

0.050 

34 

7.5 

0.133 

3.0C 

0.126 

0.130 

35 

13.0 

0.077 

3.00 

0.264 

0.260 

36 

23.0 

0.043 

3.00 

0.611 

0.590 

37 

30.0 

0.033 

3.00 

0.990 

0.960 
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TABLE  VI  -  COMPARISON  OF  OBSERVED  AND  CABUOY  CALCULATED  MOTION  OF  A 
CYLINDRICAL  WEIGHT  <L/D  -  0.75)  ON  THE  30-FEET  LENGTH  OF  BUNGEE. 

OBSERVED  CALCULATED  OBSERVED 


RUN 

PERIOD 

(sec) 

FREQUENCY 

(Hz) 

INPUT 

AMPLITUDE 

P-P 

(ft) 

OUTPUT 

AMPLITUDE 

P-P 

(ft) 

OUTPUT 

AMPLITUDE 

P-P 

(ft) 

MAGNIFICATION 

FACTOR 

(yout/yin) 

1 

30 

0.033 

6.0 

7.00 

7.36 

1.167 

2 

15 

0.067 

6.0 

5.06 

5.96 

0.843 

3 

8 

0.125 

6.0 

2.42 

2.42 

0.403 

4 

5 

0.200 

6.0 

1.07 

0.97 

0.178 

5 

3.5 

0.296 

6.0 

0.50 

0.55 

0.083 

20 

4 

0.250 

3.0 

0.25 

0.22 

0.083 

21 

5 

0.200 

3.0 

0.62 

0.61 

0.207 

22 

8 

0.125 

3.0 

1.37 

1.46 

0.457 

23 

12.5 

0.080 

3.0 

2.50 

3.22 

0.833 

24 

16 

0.0625 

3.0 

3.25 

4.22 

1.083 

25 

23 

0.435 

3.0 

3.75 

4.25 

1.250 

26 

30 

0.033 

3.0 

3.75 

3.80 

1.250 
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CABUOY 

MODEL  DEFINITION 


Figure  1.  CABUOY  (Model  Definition). 


Figure  2.  Linearization  Of  Bungee  Stress-Strain  Curve. 
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Figure  5.  Vertical  Motion  Experimental  Arrangement  In  The  100-Foot  NAVSURFWPNCEN 

Undersea  Weapons  Tank. 


Tension. 
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Figure  8.  Static  And  Dynamic  Spring  Constants  For  Multistrand  8ungee  Under  Load 


JG  CONSTANT  (LBS/FT)  FOR  30  FT  LENGTH 


4.0  • 
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Figure  9.  Comparison  Of  Dynamic  And  Static  Equations  To  Characterize  The  Multistrand  Bungee. 


OBSERVED  DATA 


Comparison  Of  CABUOY  Results  With  Observed  Measurements  For  The  Disk  System  |30-Feet  Compliance, 

6-Feet  p-p  Input  Amplitude). 


FREQUENCY  (Hz) 

Of  Obwrved  And  CABUOY  Results  For  The  Disk  System  (30-Feet  Compliance)  With  3-Feet  And 
6-Feet  p-p  Input  Amplitudes. 


Comparison  Of  Calculated  Data  On  Disk  System 


TWO  DISKS 


FREQUENCY  (Hz) 

Comparison  Of  Single  Disk  And  Two  Disk  Systems  With  3-Feet  p-p  Input  Amplitude  (30-Feet  Bungee). 


29 
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Figurt  15.  Hydrodynamic  Inertia  Coefficient  Versus  Length-To-Diameter  Ratio  For  Osciilating 

Cylinders. 

30 


I 


Comparison  Of  Cylinder  (C/d  =  0.75)  Oscillation  On  30-Feet  Bungee  With  3  And  6-Feet  p-p  Input  Amplitudes. 


Hydrodynamic  Mass  Coefficient  For  Low  Aspect  Ratio  Cylinders  For  The  CABUOY  Program  (Preliminary). 
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